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Abstract 

The early universe provides a unique laboratory for probing the frontiers of particle physics in general 
and neutrino physics in particular. The primordial abundances of the relic nuclei produced during the 
first few minutes of the evolution of the Universe depend on the electron neutrinos through the charged- 
current weak interactions among neutrons and protons (and electrons and positrons and neutrinos), and 
on all flavors of neutrinos through their contributions to the total energy density which regulates the 
universal expansion rate. The latter contribution also plays a role in determining the spectrum of the 
temperature fluctuations imprinted on the Cosmic Background Radiation (CBR) some 400 thousand 
years later. Using deuterium as a baryometer and helium-4 as a chronometer, the predictions of BBN 
and the CBR are compared to observations. The successes of, as well as challenges to the standard 
models of particle physics and cosmology are identified. While systematic uncertainties may be the 
source of some of the current tensions, it could be that the data are pointing the way to new physics. 
In particular, BBN and the CBR are used to address the questions of whether or not the relic neutrinos 
were fully populated in the early universe and, to limit the magnitude of any lepton asymmetry which 
may be concealed in the neutrinos. 
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1. Introduction 



During its early evolution the universe was hot and dense, passing brief epochs as a universal particle 
accelerator and as a cosmic nuclear reactor. As a consequence, through its evolution the entire universe 
provides a valuable alternative to terrestrial accelerators and reactors as probes of fundamental physics at 
the highest energies and densities. Several decades of progress have validated this Particle Astrophysics and 
Particle Cosmology approach to testing and constraining models of High Energy Physics and Cosmology; 
for early work see, e.g., This strategy has proven especially useful in connection with the physics 

of neutrinos (e.g., masses, mixing, number of flavors, etc.). 

Neutrinos play two different, but equally important roles in Big Bang Nucleosynthesis (BBN). On the 
one hand, electron-type neutrinos (and antineutrinos), through their charged current, weak interactions help 
to regulate the neutron-proton ratio, which plays a key role in determining the abundance of helium ( 4 He) 
emerging from BBN when the universe is ~ 20 minutes old. For example, since the He yield is largely fixed 
by the supply of neutrons available at BBN, an asymmetry between v e and D e (neutrino "degeneracy" ) will 
drive the relative abundance of neutrons up or down, thereby increasing or decreasing the relic abundance 
of 4 He. On the other hand, all flavors of neutrinos were relativistic at BBN (-few MeV > T > 30 keV), 
contributing significantly to the total density, which determines the expansion rate of the universe at that 
epoch. The competition between the universal expansion rate (the Hubble parameter, H) and the nuclear 
and weak interaction rates is key to regulating the relic abundances of the light nuclides (D, 3 He, 4 He, 7 Li) 
synthesized during BBN. 

This latter effect of (light, relativistic) neutrinos on the expansion rate also plays a role some 400 kyr later, 
at "recombination" (protons and electrons combine to form neutral hydrogen) when the Cosmic Background 
Radiation (CBR) is set free from the tyranny of electron scattering to propagate throughout the Universe. 
By influencing the age of the Universe and the size of the sound horizon at recombination, the neutrinos 
help to fix the scales of the CBR temperature anisotropies observed by WMAP and other CBR detectors; 
see, e.g., 0] and references therein. Here, however, neutrino degeneracy plays no role except, perhaps, 
by increasing the neutrino energy density and, thereby, affecting the expansion rate. This latter effect is, 
generally, subdominant. 

Since neutrinos influence the early evolution of the Universe at these two, widely separated epochs (<~ 20 
minutes and ~ 400 kyr later), the relics from BBN (light nuclides) and the temperature anisotropies imprinted 
on the CBR provide two, largely independent windows on neutrino physics. These connections and what 
we have learned from them are reviewed here. For further details and references, see |S1|S]; this review is 
largely based on these two papers. After introducing some notation in the next section, the constraints from 
the CBR are reviewed in §3. §4 provides an overview of BBN and of the current status of the comparisons 
between the observational data and the predictions of the standard model (SBBN) as well as of some general 
extensions of the the standard model (non-standard BBN). In §5 the constraints from the CBR and from 
BBN are combined to identify the allowed ranges of the baryon and neutrino parameters. We conclude in §6 
with a summary and with an identification of the successes of the standard models of particle physics and 
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cosmology and of some of the challenges confronting them. 



2. Notation 

To set the scene for the discussion to follow, it is useful to first introduce some notation. We are interested 
in three, key parameters: the baryon density, the number of "equivalent" neutrinos, and a measure of a 
neutrino-antineutrino asymmetry. 

As the universe expands, the baryon density decreases. A dimensionless measure of the baryon density is 
provided by the ratio of baryons to photons (in the CBR) . Following e ± annihilation, this ratio is preserved 
during the subsequent evolution of the universe. The parameter r\ is defined by the present (i.e., post-BBN, 
post-recombination) value of this ratio: 77 = (hb/w 7 )o; ??io = lO 10 ??. An equivalent measure of the baryon 
density is provided by the baryon density parameter, 17b, the ratio (at present) of the baryon density to the 
critical density. In terms of the present value of the Hubble parameter, Hq = lOOh km s _1 Mpc" 1 , these 
two measures are related by 

Tjio = 10 10 (n B /n 7 )o = 274n B h 2 . (1) 

In the standard model of particle physics there are three families of light neutrinos (N„ = 3) which, in 
the standard model of cosmology, are relativistic at BBN and also at recombination. The early evolution 
of the universe is "radiation dominated", i.e., the energy density is dominated by the contributions from 
relativistic particles, including the neutrinos. The universal expansion rate, as measured by the Hubble 
parameter, depends on the density: H oa p 1 ^ 2 . Any additional (non-standard) contributions to the energy 
density (such as, e.g., from additional flavors of neutrinos) will result in a speed-up of the expansion rate, 

S ee H'/H = {p'/pf 2 > 1. (2) 

Any non-standard contribution to the density may be written in terms of what would be the energy density 
due to an equivalent number of "extra" neutrinos AN,, (N„ = 3 + AN V ). Prior to e^ annihilation, this may 
be written as 

(3) 

p 43 

Thus, either S, the expansion rate factor or, AN„, the number of equivalent neutrinos, provide equally good 
measures of the early universe expansion rate. While it is easy to imagine extra contributions to the energy 
density from new physics beyond the standard model, it must be noted that it is possible for AA^ to be 
negative, leading to a slower than standard, early universe expansion rate (S < 1). For example, models 
where the decay of a massive particle, produced earlier in the evolution of the universe, reheats the universe 
to a temperature which is not high enough to (re)populate a thermal spectrum of the standard neutrinos 
(Trh < 7 MeV), will result in AJV„ < and S < 1 0. 

For any neutrino flavor i, an asymmetry ("neutrino degeneracy") between the numbers of vi and T>i, 
relative to the number of CBR photons, can be quantified by the net lepton number L,;, the neutrino 
chemical potential pi or, by the dimensionless degeneracy parameter £j = Hi/T: 
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Although we are interested in lepton asymmetries which are orders of magnitude larger than the baryon 
asymmetry (B ~ rj < 10~ 9 ), the values of £j [i = e,/j,,r) considered here are sufficiently small 0.1) so 

that the "extra" energy density contributed by such degenerate neutrinos is negligible. 

AN„&) = ^ (%» + ^(&)*<0.01. (5) 

7 7T 7 7T 

In this case, the results to be presented below for £ ^ will correspond to = 3 (S = 1). In fact, if the three 
active neutrinos (u e , , v T ) mix only with each other, all individual neutrino degeneracies will equilibrate 
via neutrino oscillations to, approximately, the electron neutrino degeneracy before BBN begins [Sj- Thus, 
the magnitude of the electron neutrino degeneracy constrained by BBN is of special interest when limiting 
the total net lepton asymmetry of the universe: L 3L e ; AN„(£) « 3AN„(£ e ). 

For the standard models of particle physics and cosmology, AiV„ = (S = 1) and £ = £ e « £^ « £ T = 0, 
and the value (range of values) of rj identified by BBN and the CBR should agree, restricting the allowed 
deviations from zero of AN„ and/or £ e . 

3. CBR 
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Figure 1: The CBR temperature anisotropy angular power spectrum for three choices of the baryon density, 
Q,^h 2 = 0.018, 0.023 (best fit), 0.028, from bottom to top respectively near I 200. The data points with 
error bars are from Spergel et al. 0]. 

In Figures 1 and 2 are shown the CBR temperature anisotropy angular power spectra for different choices 
of the baryon density (Fig. 1) and of N„ (Fig. 2). Non-zero values of AiV„ change the energy density in 
radiation, which shifts the redshift of the epoch of equal matter (Cold Dark Matter and Baryons) and 
radiation densities. This results in changes to the angular scales and the amplitudes of the "acoustic" peaks 
in Figures 1 & 2; see, e.g., ,SJ and further references therein. WMAP is a much more sensitive baryometer 
than it is a chronometer. While the best fit values for the baryon density and N„ are r]io = 6.3 and N^ = 
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2.75 (S = 0.98) respectively, the 2a range for the baryon density parameter is 5.6 < r/ 10 < 7.3, whereas for 
N„ it is 0.9 < N„ < 8.3 (0.81 < S < 1.36) 0. Thus, although the WMAP best fit value of N„ is less than 
the standard model value of 3, it is clear that this difference is not at all statistically significant. It will 
be interesting to see if the much tighter CBR constraint on the baryon density parameter (~ 6 — 8 %) is 
consistent with the value of this parameter identified by SBBN. 
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Figure 2: The CBR temperature anisotropy angular power spectrum as in Fig.^ now for four choices of the 
equivalent number of neutrinos N„ = 1, 2.75 (best fit), 5, 7, from bottom to top respectively near I ss 200. 
The data points with error bars are from Spergel et al. 



4. BBN 

Since the relic abundances of D, 3 He, and 7 Li produced during BBN are rate limited (nuclear reaction 
rates), each of these nuclides is a candidate baryometer. Among these, Deuterium is the baryometer of 
choice for several reasons. The BBN-predicted abundance of D is sensitive to the baryon density (j/d = 
10 5 (D/H)p cx rj^ 1 - 6 ). The post-BBN evolution of D is simple in that whenever gas is incorporated into stars, 
D is completely destroyed. Thus, observations of the deuterium abundance anywhere, at any time, provide 
a lower bound to the primordial D abundance (and, therefore, an upper bound to rf) . It is expected that if 
D can be observed in regions which have experienced minimal stellar processing, the deuterium abundance 
inferred from such data should be very close to the BBN abundance. 

The good news is that there are data from observations of neutral D and neutral H in high redshift, 
low heavy element ( "metallicity" ) abundance QSO absorption line systems (QSOALS); see Figure [21 As 
may be seen from Fig. |3 the bad news is that there are only five such systems with good enough data to 
derive meaningful D abundances 9 . And, even for these, specially selected targets, there is the possibility 
of confusion between D and H absorption spectra which are identical, save for the wavelength/velocity shift 
between them. That is, small amounts of hydrogen at the "wrong" redshift" (interlopers) can masquerade as 
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deuterium. Further, since the hydrogen absorption in such systems is saturated, it is often difficult to identify 
the number of systems which contribute to the total absorption and this can lead to errors in the inferred 
amount of H in determining the D/H ratio. With these caveats in mind, it is important to understand 
that systematic, rather than statistical uncertainties may dominate the error budget. Indeed, for the data 
summarized by Kirkman et al. JOJ and shown in Fig. x 2 exceeds 16 for 4 degrees of freedom! Following 
Kirkman et ai, the error bars are inflated here to account for this excessive dispersion, and a primordial 
abundance j/d = 2.6 ±0.4 is adopted. For SBBN this estimate for the abundance of primordial D corresponds 
to f]io = 6.1+q 5 . This is in excellent (spectacular!) agreement with the completely independent estimate 
above in §3 from the CBR. In Figure 01 are shown the likelihood distributions for r\ derived from the CBR 
(WMAP) and from SBBN (D). 

3 He is a much less useful baryometer than is D. In the first place, its BBN-predicted abundance is less 
sensitive to the baryon density parameter ( 3 He/H oc ?7~ ' 6 ). In addition, as gas is incorporated into stars and 
the stellar-processed material is returned to the interstellar medium when they die, 3 He is produced, destroyed 
and, some survives. This complicated history makes it much more difficult to account for the post-BBN 
evolution of 3 He. Finally, 3 He is only observed within the Galaxy ^H] and, therefore, its abundance samples 
only a limited range in metallicity. Nonetheless, while there is a clear oxygen abundance gradient with 
location within the Galaxy (higher in the center, lower in the suburbs, indicating more stellar processing in 
the interior), the 3 He abundance shows no such gradient (either with position or with metallicity). However, 
the Bania, Rood, and Balser ^H] recommended value of j/3 = 10 5 ( 3 He/H) = 1.1 ± 0.2 corresponds to 
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Figure 3: The deuterium abundance (by number relative to hydrogen), = 10 5 (D/H), derived from high 
redshift, low metallicity QSOALS [§] (filled circles). The metallicity is on a log scale relative to solar; 
depending on the line-of-sight, X may be oxygen or silicon. Also shown is the solar system abundance (filled 
triangle) and that from observations of the local ISM (filled square). 
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T] W « 5.6^17, in broad agreement with the SBBN-D and CBR-WMAP determinations. So far, so good. 

In contrast to deuterium (and to 3 He), the abundance of 4 He has increased in the post-BBN universe 
as stars have burned hydrogen to helium (and beyond). Therefore, to avoid model-dependent evolutionary 
uncertainties, it is best to concentrate on determining the 4 He abundance (mass fraction, Yp) in the most 
metal-poor sample available and to let the data speak for themselves concerning any correlation between the 
helium and heavy element abundances. The best such data come from observing the emission lines formed 
when ionized helium and hydrogen capture electrons in regions of hot, ionized gas (H n regions). There exists 
at present, thanks largely to the work of Izotov and Thuan a very large sample of helium abundance 
determinations in low metallicity, extragalactic H n regions. This newer, more uniform data set complements 
earlier, more heterogeneous samples ^21- For the WMAP estimate of the baryon density, including its 
uncertainty, the SBBN-predicted primordial abundance of 4 He is Yp = 0.2482 ± 0.0007. Unfortunately, none 
of the Yp estimates IT2*] agree with the SBBN prediction, all being low by roughly 2a. Indeed, from 
their 2004 sample of 82 data points Izotov & Thuan [II derive such a small uncertainty, that their central 
value is low by nearly 6cr! For the numerical results presented in §5 below, a primordial mass fraction Yp = 
0.238±0.005 Q3] is adopted. 

The 4 He abundance determinations are most likely examples of extremely precise, yet inaccurate, de- 
terminations of an important cosmological parameter. It has long been known that there are a variety of 
systematic uncertainties which are likely to interfere with an accurate Yp determination. In a very recent, 
detailed study of some of these identified systematic uncertainties, Olive & Skillman |14j suggest the true 
errors likely exceed previous estimates by factors of 2-3 or more (AYp as 0.013). Given such large uncer- 




7 7lO 



Figure 4: The likelihood distributions for the baryon density parameter 7710 inferred from the CBR (WMAP) 
and from SBBN and deuterium; see the text. The broader distribution, centered at 7710 ~ 6.1 is from SBBN 
and the narrower distribution, centered at 7710 ~ 6.3 is from the CBR. 
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taintics, it is not surprising that the extant 4 He abundance data are, within the errors, consistent with the 
predictions of SBBN and the CBR (and/or SBBN plus D) determined baryon density. Nonetheless, it might 
be premature to ignore this challenge to the standard models of particle physics and of cosmology. Perhaps 
the tension between D and 4 He is an early warning of new physics beyond the standard models. Before 
pursuing this option in the next section, 7 Li is considered here. 

As with 4 He, there is conflict in the comparisons between the SBBN predictions and the 7 Li relic abun- 
dance estimates derived from observations. Here, too, the potential for systematic errors looms large. 7 Li is 
produced in the Galaxy by cosmic ray spallation/fusion reactions and observations of super-lithium rich red 
giants provide evidence that (at least some) stars can be net producers of lithium. Therefore, to infer the 
BBN yield of 7 Li, the data should be limited to those from the most metal-poor halo stars in the Galaxy. 
For the WMAP baryon density, the SBBN expected 7 Li abundance is [Li] P = 12+log(Li/H) P = 2.65±o °g . 
In contrast, for a selected data set of the lowest metallicity halo stars, Ryan et al. |I5| derive a primordial 
abundance of [Li]p k, 2.0 — 2.1. In deriving the stellar lithium abundances, the adopted stellar temperature 
plays a key role. When using the infrared flux method effective temperatures, studies of halo and Galactic 
Globular Cluster stars JHj suggest a higher abundance: [Li]p = 2.24 ± 0.01. Very recently, Melendez & 
Ramirez |17j have reanalyzed 62 halo stars using an improved infrared flux method effective temperature 
scale, confirming the higher lithium abundance; they find [Li]p = 2.37±0.05. If this were the true primordial 
7 Li abundance, then the SBBN value of the baryon density parameter would be 7710 = 4.5 ± 0.4, in conflict 
with the CBR- WMAP and/or SBBN-D estimates. Indeed, all of the current observational estimates of the 
abundance of primordial lithium arc significantly lower than the SBBN expectation. 

As with 4 He, the problem may be traced to the astrophysics rather than to the cosmology. Since the low 
metallicity halo stars used to estimate the primordial abundance of lithium are the oldest stars in the Galaxy, 
they have had the most time to modify (by dilution and/or destruction) their surface abundances. While 
mixing of the stellar surface material with the interior would destroy or dilute prestellar lithium, the very 
small dispersion among the observed values of [Li] derived from the lowest metallicity halo stars suggests this 
effect may not be large enough to bridge the ~ 0.3 dex gap between the observed and CBR/SBBN-predicted 
abundances; see, e.g., |18j and further references therein. 

5. CBR And BBN Combined 

In contrast to D (and 3 He and 7 Li), 4 He is an insensitive baryometer, but its primordial abundance is a 
useful, early universe chronometer. If the standard model expansion rate is changed (5^1, AN„ ^ 0), 
this will affect the neutron-proton ratio at BBN and change the SBBN-expected value of Yp. The current 
conflict between the SBBN-predicted and the observationally inferred values of Yp requires a slowdown in 
the early universe expansion rate (S < 1; AN V < 0). A joint BBN fit to the observationally inferred D and 
4 He abundances suggests that 7710 ~ 5.7 and S « 0.94 (AiV„ ks —0.70) can relieve the SBBN tension between 
D and 4 He. However, it can be seen in Figure |S] that while these values are entirely consistent with the 
constraints from the CBR, N„ = 3 is consistent with both BBN and the CBR at ~ 2a 0. Nonetheless, this 
combination of parameters does not resolve the conflict with 7 Li. Although a slowdown in the expansion 
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rate has the effect of increasing 7 Li (more time for production), this is compensated by the somewhat lower 
baryon density (slower reaction rates), which has the opposite effect. The result is that for the choices 
of S and ij which resolve the conflicts between D and 4 He (and between WMAP and 4 He), the predicted 
primordial abundance of 7 Li is [Li]p w 2.62 ± 0.10, still some 0.2-0.3 dex higher than that inferred from the 
data. 
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Figure 5: The lcr, 2a, and 3a contours in the AiV„ - 7710 plane (from Barger et al. [5]) consistent with the 
WMAP CBR data and the BBN predicted and observed abundances of D and 4 He. The cross marks the 
best fit point; see the text. Note that the upper horizontal axis is for wb = f^e/i 2 . 

Another example of new physics with the potential to resolve the D - 4 He conflict while leaving the CBR- 
D agreement unaffected is a neutrino- antineutrino asymmetry (£ = £ e ^ 0); see [H] and references therein. 
Through its effect on the neutron-proton ratio at BBN, such an asymmetry can change the BBN-predicted 
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Figure 6: The la and 2a contours in the £ e — r/io plane (from Barger et al. 0) consistent with the WMAP 
CBR data and the BBN predicted and observed abundances of D and 4 He. The cross marks the best fit 
point; see the text. 
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4 He abundance. For £ e > 0, there are more v e than v e and this drives the n-p ratio down, decreasing the BBN- 
predicted 4 He abundance. For £ e < 0.1, the extra energy density contributed by these degenerate neutrinos 
is small, and N„ = 3 remains a good approximation. As a result, such a lepton asymmetry has negligible 
effect on the CBR temperature anisotropics and the good agreement with the WMAP data is unaffected. 
The effects of such an asymmetry on the predicted abundances of D, 3 He, and Li are subdominant to that 
on 4 He. The best fit parameter choices which resolve the tension between D and 4 He, while preserving the 
good agreement with the WMAP data are 7710 = 6.2 (Q,B,h 2 = 0.23) and £ e = 0.044. In Figure are shown 
the la and 2a contours consistent with the WMAP data and with BBN (D and 4 He). 

Finally, we note that if all three parameters (77, 5, and £) are allowed to be free, much larger ranges in 
them will remain consistent with BBN (D and 4 He) and with the CBR; see, e.g., [H] and the more recent paper 
of Kneller and Steigman |19| . Fixing the D and 4 He abundances, Kneller & Steigman find two approximate, 
but quite accurate, BBN relations among these three parameters. 

590(5 - 1) « 1161710 - 697, (6) 

and 

145£ e w 106(5 - 1) + 6.31. (7) 

Consistent with the WMAP CBR data, values of AN„ and £ e in the ranges —2 < AN,, < + 5 and 
—0.1 i$ £ e ~ + 0.3 are permitted 0. However, even with this freedom it is still not possible to recon- 
cile the BBN-predicted and the observed relic lithium abundances [FQ. For the values of these parameters 
which are consistent with BBN and the CBR, [Li]p ~ 2.6. 

6. Summary and Conclusions 

BBN and the CBR probe the evolution of the Universe (and its constituents) at two, widely separated epochs 
in its early evolution. Confronting the predictions of BBN and the CBR with the relic abundance and WMAP 
data enables independent tests of the standard models of particle physics and cosmology. Qualitatively, the 
standard models pass these tests with flying colors, permitting BBN and CBR constraints to be put on 
new neutrinos physics (N„ 7^ 3?; £ e 7^ 0?). When considered in quantitative detail however, there are some 
challenges to the standard models at the ~ 2a level. Many would take this as evidence for success and 
declare victory. However, if these tensions are taken at face value, they might be alerting us to problems 
with the astronomy, the astrophysics, the cosmology, the particle physics or, combinations of them. It is 
not unlikely that the apparent conflicts may result from the data, its analysis, and/or the extrapolations to 
the early universe. While the community awaits the new surprises to be encountered at the LHC, a Linear 
Collider, or the next generation of terrestrial or space-based telescopes, it should be kept in mind that these 
challenges could be pointing the way to new physics, especially new neutrino physics, beyond the standard 
models of particle physics and/or cosmology. 
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